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An analytical study is made of the potential distribution 
about an electrode on a stratified earth, from the point of 
view of developing special formulas appropriate for 
computation at small and large distances from the elec- 
trode. Power series in the distance for small distances and 
simple asymptotic expansions for large separations are 
derived first for a two layer earth, from integral expressions 
for the potentials. For graphical integration a transfor- 
mation is derived which makes the integrals converge 
very rapidly. For a many layer earth a general method 
for deriving the power series for small distances is indi- 
cated. For large distances, general expansions are given in 
closed (Eq. (23)) and asymptotic forms (Eq. (24)) applying 
to all types of conductivity variation. Explicit expressions 


for the power series are given for a three layer earth with 
the middle layer of thickness equal to and twice that of 
the upper layer. Closed and very rapidly converging series 
of Hankel function expansions are then developed for the 
third layer having zero or infinite conductivities, with the 
thickness of the middle layer being equal to, twice and 
three times that of the upper layer. Thirty numerical 
cases have been computed and plotted as curves to 
illustrate the fundamental formulae developed. In an 
appendix general proofs are given for the asymptotic 
expansions of not only the type of integrals occurring in 
the present problem but also for integrals irivolving Weber 
and Hankel functions as well as the Bessel functions, and 
where these functions are of non-zero orders. 


INTRODUCTION 


GREAT deal of work has been done on 

the problem of the potential distribution 
about a point electrode placed on the earth’s 
surface, the earth being supposed to consist of 
two or more plane faced layers of different 
conductivities. Thus as long ago as 1873 Max- 
well! set up the system of images with which to 
treat the problem. Recently the problem has 
again attracted attention. Curves for the po- 
tential distribution, but without the analysis, 
were given by Weaver.® Following these, formu- 
lae in terms of images, but without explicit deri- 
vation, were published by Hummel.’ Stefanesco 
and Schlumberger* then gave a derivation of 


1J. C. Maxwell, Electricity and Magnetism, 1873. 

2 W. Weaver, Amer. Inst. Min. and Met. Eng. Tech. 
Pub. 121, 1928. 

3 J. N. Hummel, Zeits. Geophysik 5, 89, 228 (1929). 

4S. Stefanesco and C. and M. Schlumberger; J. de 
Phys. et Rad. 1, 132 (1930); cf. also L. J. Peters and 
]. Bardeen; Physics 2, 103 (1932). 


Hummel’s formulae by a direct solution of 
Laplace’s equation in terms of Bessel functions. 
Finally, Roman® published extensive numerical 
tables for the computation of the potential 
distribution, for a two layer earth. 

Although these investigations have provided 
a number of numerical solutions nothing seems 
to have been given in the sense of a detailed 
analytical discussion of the formal solutions. 
In connection with the problem of the flow of 
underground waters in artesian wells,® the writer 
had occasion to study in detail the potential 
distribution about an electrode in an insulated 
place faced conducting slab. Because of the great 
geophysical interest in the more general problem 
of several conducting slabs it was thought 
worthwhile to generalize the methods of the 
previous analysis and apply them to the case of 
several conducting slabs. In the present paper 
we shall first present the generalizations applying 


5T. Roman; U.S. Bur. Mines, Tech. Paper 502, 1931. 
®M. Muskat, Physics 2, 329, 1932. 
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to the simpler problem of a plane faced surface 
layer overlying a semi-infinite bed of conductivity 
different from that of the surface layer and then 
extend them to the case of a three layer earth. 
The essential results will be the expressions of 
the potential distributions near the electrode as 
power series in the radial distance from the 
vertical axis through the electrode and as 
asymptotic expansions in the reciprocal of this 
distance for large distances from the electrode 
axis. 


I. AN ELECTRODE IN A Two LAYER EARTH 


We shall base our fundamental formulae upon 
the solution of Laplace’s equation in terms of 
Bessel functions. Our coordinate system will 
be chosen as indicated in Fig. 1. The problem 
then will be to find solutions ©, and #2 of 


vb =0 ) 
such that 
:2=a/2; 
(1) 
OP, OP, OP, a 
—=0 :2=0; 
Oz 0z 2) 


and so that , have a simple pole at (0, /f). 
Before proceeding with the solutions it is con- 
venient to take the unit of length equal to twice 
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the thickness of the surface layer, and introduce 
the dimensionless variables: 


w=z2/a; p=r/a; x=f/a. (1a) 


With regard to the conductivity constants 
a1, 02 it will be useful to introduce both the 
symbol k, where: 


k 


and the quantities y defined by: 


y=ctnh (2/01) : o2/o,>1: 


=tanh-" (a2/01) : o2/0,<1 (1c) 
so that: 


=—e 7; :k=e7*7 : o2/0,<1. (1d) 


Although the final results suitable for numer- 
ical discussion may be most conveniently ex- 
pressed in terms of k, the derivation of the 
formulae of interest may be more directly 
developed by the use of the symbol y, and will 
therefore be employed in most of the formal 
steps of the analysis, even though this will 
necessitate the separate treatment of the cases 
o2/0;>1 and o2/0,<1. 

Since Stefanesco and Schlumberger‘ have 
already given the method of solving Eq. (1) by 
means of Bessel functions we shall state the 
results only. They are: 


© J((2pa) cosh 2aw sinh (a+y—2ax)da ) 


o;/o2<1 


0 


cosh 


© J,(2pa) cosh 2ax sinh (a+y—2aw)da 
2f 
0 


w=}=x sinh vf 
0 


2 
cosh 


Jo(2pa) cosh 


cosh (a+7) 


Jo(2pa) cosh 2aw cosh (a+y—2ax)da ) 
w=x=} 
0 


sinh (a+~7) 


J((2pa) cosh 2ax cosh (a+~y—2aw)da 
:%,=2 


0 


w=4=x cosh y 
0 


. f (3) 
sinh (a+~7) 


J(2pa) cosh 2axe~@vVada 


sinh (a+7) 
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Fic. 1. Diagrammatic representation of an electrode, at 
(0, f) in a two layer earth. 


Eqs. (2) and (3) are essentially equivalent to 
those of Stefanesco and Schlumberger, being 
only somewhat more symmetrical than the 
corresponding formulae already in the literature. 
The symmetry has been attained by the intro- 
duction of the quantity y rather than the k of 
Eq. (ib). On the other hand, the analysis of 
Stefanesco and Schlumberger really went no 
further than the demonstration that their equiva- 
lents of Eqs. (2) and (3) could be evaluated in 
the form of the series of images already given by 
Hummel. Here, however, we shall begin with 
these equations and derive from them separate 
expansions suitable for numerical treatment 
when p is small and when » is large. 


A. Expansions for small values of p 


To study these solutions for small value of p 
we shall expand the denominators in Eqs. (2) 
and (3) as: 


sech (a+7) = 23: n( —1)"e~ 
(4) 


esch = 23 ne- ent 
0 


We might now integrate term by term on 
using the relation: 


= (a+ (4a) 


obtaining at once Hummel’s series of images 
with terms of the form 


These could then be individually expanded as 
power series in p and summed. It will be more 
convenient, however, to expand Jo(2pa) and 
then integrate the double series. For purposes 
of brevity we shall restrict ourselves to the 
problem of most practical interest—the potential 
distribution on the earth’s surface due to an 
electrode on the surface, i.e., the case where 
w=x=0. Furthermore, since we shall consider 
only the function ®,, we shall drop the subscript. 
Carrying out, then, the procedure indicated, and 
returning to the notation of Eq. (1b) we find: 


If k=0, i.e., if there is no conductivity contrast 
between the two media, then = 1/p, as it should 
be for the surface potential due to an electrode 
embedded in the surface of a semi-infinite 
homogeneous medium. If, however, the con- 
ductivity contrast is nonvanishing, |k| >0, the 


series have to be summed. 

In doing this it is convenient to take out and 
sum separately the terms for m=0 (the constant 
terms) and those for »=1 and 2, corresponding 
to the images nearest the electrode. Eq. (5) can 
then be rewritten as: 


b= 1/p—2 log (1—k) (1+ (4+ (6)' 
1 


7 The constants independent of p have to be retained in order to make the above ® connect with that to be derived 
later for large p from the same integrals of Eq. (2) and Eq. (3). 
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where 


Cm = fn(k) = Som (7) 
3 


The series for f,,(k) evidently converges very 
rapidly even if k~1. Its value, in fact, will 
always be of the order of k*/3?"*'. Hence, since 
Cm=O(1/m)' for m large, the radius of con- 
vergence of the series in Eq. (7) is given by 
|p| =3.8 Thus, whereas the direct image series 
converges essentially as k"/n, the expansion of 
Eq. (6) converges in p as p*"/3™. 

It is of interest to note that if only the terms 
for m=0 are separated from Eq. (5) it will 
reduce, for the case o2.=0, (k=1) to: 


&=1/p—2 log (1—k) +25-m(—1) (2m-+1)p 
1 


where ¢ is the zeta-function. This is equivalent 
to Eq. (7) in the analysis of the problem of the 
insulated conducting slab previously treated by 
the author,® as of course it should for o2.=0. 

Before illustrating the use of Eq. (6) numeri- 
cally, we shall develop the corresponding ex- 
pressions suitable for large values of p. 


B. Asymptotic expansions for ® for large p 

To study our potential function for large 
values of p it is convenient to return again to 
the integral expressions of Eqs. (2) and (3). 
However, we cannot simply expand sech (a+ 7) 


and csch (a+) as in Eq. (4) and then apply the 
asymptotic expansion for Jo(2pa). For as the 
asymptotic expansion of Jo(2pa) involves terms 
as 1/(a)" with m>1 the resulting integrals 
would diverge at a=0. Neither could we ex- 
pand directly the image series of terms like 
e-2"7/(p?-+-n2)!, since even if p is large, there will 
always be terms for which 1> p. 

Before, however, proceeding with the deriva- 
tion of a real asymptotic expansion for ® it 
may be of interest to note that if one wishes to 
integrate Eqs. (2) and (3) graphically a more 
suitable form, especially for large p, may be 
obtained as follows. First, Jo may be split into 
the two Hankel integrals as: 


2J0(2 pa) = (2 pa) +119 (2 pa).® 


The two separate integrals formed in this way 
(from the original integrals of Eqs. (2) or (3)) 
may then be taken with contours in the complex 
a plane; the integral with H,"?(2pa) may be 
carried over the infinite first quadrant arc and 
down the + imaginary axis, and the integral 
with Hy (2pa) may be carried over the infinite 
4th quadrant arc and up the — imaginary axis. 
Noting then that since y>0 there are no poles 
in these contours, and using the relation: 


(iz) = — ( —iz) = (—2i/r) Ko(2), 
Eqs. (2) and (3), in the notation of Eq. (1b) reduce to: 


2 Ko(pu) cos wy cos xp du ) 
0 1+k?—2k cos 


cos xu(cos cos u(1—w))du 


(8) 


2 
=1:6=-(1+k 
we 


1+k?—2k cos 


Since Kg vanishes exponentially for large arguments Eq. (8) is convenient for graphical 
integration for large p, unless k~+1(o2~0). The case k~—1, however, will cause but little 


® This radius of convergence could be advanced by 
separating out still further image terms, as those corre- 
sponding to n=3, and m=4. On the other hand, without 
any separation, the series of Eq. (5) has a radius of con- 
vergence of |p| =1. It may be noted, too, that Roman’s® 
rather lengthy and complex analysis is really equivalent 


to separating out the terms for »=1---5 and expressing 


the remainder as a series in (p/10)? with numerical coeffi- 
cients. 

® This, and the other relations involving the Bessel or 
related functions which will be used below, may be found 
in Theory of Bessel Functions by Watson. 
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difficulty since the minimum of the denominators will occur first at 4= a and even if p=1, 
Ko(pu)/(1+)? <3 for k as small as —0.90. 

On the other hand, Eq. (8) is unsuited for the derivation of an asymptotic expansion. For if 
Ko(pu) be expanded under the integral sign the same difficulties will arise as in the case of the 
expansion of Jo in the original Eqs. (2) and (3). We have therefore to resort to quite a different 
method; and since we shall be able to express %, in a rigorous closed series form (cf. Eqs. (11) 
(12) below) we shall treat the more general case where w, x0, but carry through the details only 
for o;/¢2<1 and drop the subscripts. Thus by the “‘partial fraction’”’ rule we have: 


cosh 2aw sinh (a+y—2ax) 


= Ln 
cosh (a+7) 


cosh 27(w+x) cos nx(w+x)+cosh 2y(w—x) cos nr(w—x)} 
+nr/2{sinh 2y(w+x) sin nx(w+x)+sinh 2y(w—x) sin nx(w—x)}}. 


(9) 


4 


We now apply the relations: 


= 


0 


and then integrate Eq. (2) with respect to a, making use of Eq. (4a), finding: 


cos mn py 
6=2> n  fcosh 2y(w+x) cos rn(w+x)+cosh 27(w—x) cos 
0 (1+y*)! 
r (10) 
sin mnpy 
+ {sinh 2y(w+x) sin rn(w+x)+sinh 2y(w— x) sin | 
0 (1+y*)! 
We now observe that: 
0 


where Ho denotes the Struve function of zero order and Yo is the Weber function, in Watson’s 
notation. We, therefore, have finally: 


oi/oz<1; k<O; 


=15-n{cosh w(2y—(2n+1)xi) cosh x(2y— (2n-+1)xi)| 


(11) 
{To(2yp—(2n+1) mip) — Yo(2yp—(2n+1)mip)} 
and similarly for o1/02>1 :k>0; }=w; 
2w(y—ni) cosh 2x(y—nz1)}- { — 2nmip) — Yo(2yp—2nzip)}. (12) 


Although these equations really contain the final result of our analysis for this section we shall 
restrict our further discussion to the special case of most practical interest: w= x=0. We have then, 
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0 


oi/o2>1:k>0; w=x=0; 
Jo(2pa) ctnh (y+a)da=2r> n{ Ho(2yp—2nmip) — Yo(2yp—2nzip)}. (12a) 


Making use now of the asymptotic expansion of the Struve function Hy, which is given by 


I*(m+}) 
= Yo(s)-+— (—1)* (13) 


m=0 
and then, summing the series with respect to , we get finally: 
o;/02<1: k<0; w=x=0; 
1-1 (—1)"l2(m+3) a 


tanh y+0(p-??-!); 14 


o102>1:k>0; w=x=0; 


=— t 
T T'(2m+1) ayn (15) 


Since tanh y and its derivatives are always less than 1 except for the numerical factors developing 
on successive differentiations, and as I?(m+})/T'(2m+1)~1/2?"(m)! for large m, the series of 
Eq. (14) will be quite suitable for numerical purposes even though in a strict sense it is a divergent 
series. Eq. (15), on the other hand, is not satisfactory when y is small (k~1, o:/¢2>>1) unless p 
is so large that even yp is large. For moderate values of p and small values of y we may derive a 
more satisfactory series by separating out in Eq. (12a) the term with n=0. Carrying through the 
rest of the analysis as before we get: 


w=x=0:0;/o2>1; 


where the values of Hy) and Y) may be taken from the tables of Watson. 

Although Eqs. (14) and (15) are rather compact analytically their physical significance may be 
somewhat clearer if expressed explicitly in terms of the conductivity ratios o;/o2. Carrying out the 
differentiations for the first four terms we find that both for o:1/¢2>1 and o;/02<1 the potential for 
x= w=0 has the form: 


3 
=— (: -=)- “\(= -2) += (1-5) (17-607 | . (16) 
o2p 4p? 16p* 64p° 


Before entering now upon a numerical discussion of Eqs. (14) and (15) we must point out that 
they are not to be used in the limiting cases where y=0. Thus, in the case of Eq. (14) where y=0 
means that o,/02=0 we find at once the fallacious result that @ vanishes as y approaches zero. 


13) 


14) 
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The difficulty, however, is only apparent. For it is to be remembered that Eqs. (14) and (15) are 
only asymptotic expansions and the fact that 6(y=0)=0 simply means that when y=0 the asymp- 
totic expansion for # has a leading term /ess than 1/p".'° And in fact this appears automatically 
when we put y=0 in the rigorously valid Eq. (11a), before we develop the asymptotic expansion, 


for we have then: 


=0) = mip) — Yo((2n+1)rip)}. (1b) 


For imaginary arguments we may apply the relation: 


H (iz) +Ho(—iz) — Yo(iz) — Yo( —iz) = (4/m)Ko(2), (17) 


where Ky is the ordinary Hankel function, so 
that Eq. (11b) becomes: 


(0;/02=7 =0) =4 Ko(ump); (Alc) 


as the Hankel functions vanish exponentially 
for large arguments we see that #(y=0) has 
indeed an exponential leading term, which from 
the point of view of an asymptotic expansion is 
vanishingly small compared to a term as 1/p”. 

In the case of Eq. (15) the situation is just 
reversed. There, when y~0, which means that 
the lower medium is a perfect insulator (o2=0), 
& becomes infinite. 

This difficulty, on the other hand, is real from 
an analytical point of view, as it is known that 
the potential due to the infinite series of images 
for this special case diverges. The problem then 
is to separate out the diverging term as an 
infinite constant and then use the finite variable 
part which does have physical significance." 
This may be readily done by returning to Eq. 
(12a) which shows that when y~0O the term 
with »=0 becomes infinite. Separating out this 
term, applying Eq. (17) and noting that H1)(0)=0, 
we find: 


=0=o2) = —2 log 2y—1.1544 


+2 log Ko(2nmp). (18) 
1 
Thus, the first term on the right contains the 


10 As will be seen later, this point will arise also in the 
general case of a many layered earth. 

1! This question arises also in our series expansion when 
e is small, for when k~1 the term —2 log (1—&) in 
Eq. (6) becomes infinite. 


diverging element of & (and in a manner equiva- 
lent to that in Eq. (6)) and the remainder gives 
the variable part of physical interest. And in 
fact, this variable part is actually identical 
with that given previously by the author (cf. 
reference 6, Eq. (10)) for the special case of an 
insulated slab. It may be noted, too, that only 
in this case—where the upper layer is insulated 
—does the potential at large distances from the 
electrode vary logarithmically with the distance. 
For finite conductivities of the lower layer and 
large p, the leading term in ® is proportional to 
1/p. The numerical rate of decrease of © with p, 
however, increases continuously from the loga- 
rithmic rate for o.=0 to that of 1/p for o.2=0; 
and finally to an exponential decrease for o.~ ~. 
Hence, as will be seen below (cf. Fig. 3) &p 
increases with p for o2/0,;<1 but decreases with 
p when o2/0,>1. 


C. Numerical discussion 


Since a number of numerical cases for the 
problem of this section are already to be found 
in the literature we shall confine our numerical 
discussion to several illustrations of the formulae 
we have derived by our analysis. For these we 
have chosen the cases where the conductivity 
ratios have the values: 


o2/o,=0.01, 0.20, 1.0, 5.0, 


Now it will be recalled from Eqs. (6) and (7) 
that for small values of p(p<3) our series 
expansions involve the constants c,, and the 
functions f,,(k). Since the c,, are independent of 
k they need be computed only once and they 
will apply to all values of k. The first 10 values 
are: 


la) 
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m: 1 2 3 4 5 6 7 8 9 10 
Cm: 0.5000 0.3750 0.3125 0.2734 0.2461 0.2256 0.2095 0.1964 0.1855 0.1762 


As to the functions f,,.() the computations are for 2 >0 developed as a Taylor series in y, after 
almost as simple. For it is evident from their which Eq. (17) was applied to the first terms of 
definition that the series for f,,(k) will converge this expansion.” For o2/¢,=0.20, Eq. (15a) was 
quite rapidly and in the most unfavorable casé found very satisfactory. o2/0,=1.0 corresponds, 
where k=.+1 the values may be deduced from of course, to a homogeneous earth and hence 
known series. In Fig. 2 we have plotted the first gives #=1/p everywhere, as also follows from 
four values of f,,.(k) for the range -—1=k=+1, Eq. (16). The case o2/0:,=5.0 was readily 
which includes all positive values of the ratio computed by Eq. (16), and that for o2/0o,=% 
o2/0;. Because of the very rapid fall in the value was treated by the very rapidly converging 
of f,, with increasing m they have been multiplied _ series of Eq. (11c). 
by successively increasing factors so that they The results of these calculations are plotted in 
could be plotted on the same scale. The very Fig. 3. To avoid the overlapping of the curves 
fact that this has been necessary shows perhaps and the consequent masking of the effects of 
most clearly the rapid convergence of the series the variation of o2/0:, the function @p has been 
in Eq. (6). plotted rather than ® itself. Since curves of this 
For large values of p(p>3), the chosen values type are already available in the literature we 
of o2/c, serve to illustrate the various asymptotic shall not discuss them further. However, we do 
expansions that have been developed for this wish to point out the practical value of the 
purpose. Thus, for o2/o,=0.01, Eq. (12a) was formulae developed here with respect to the 
used with the term »=0 computed separately —~_ The result of these steps was to give, first, Eq. (18) 
from the tables of Watson, and with the terms and then a small correction term proportional to y. 
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k 
Fic. 2. Plot of functions f,,(k), defined by Eq. (7). 
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Fic. 3. Surface potential distribution in a two layer earth 
o,=conductivity of surface layer, o2=conductivity of 
lower layer. p = (radial distance from electrode) /(twice the 
thickness of upper 


layer). 

calculation of such curves. For in the compu- 
tations by Eq. (6) one may drop the series after 
the first four terms and introduce an error of 
less than 0.6 percent in @p even with p= 2.0 and 
|k| ~1, for o2/0,<1. For o2/0,;>1 the errors are 
somewhat larger, but if necessary a few more 
terms may be used, each additional term re- 
ducing the error at least by the factor 4/9. On 
the other hand, as the potential distributions 
must be used graphically for practical applica- 
tions very high accuracy will seldom be required. 
Finally, it may be noted that the series of Eq. 
(6) is an alternating series and hence permits a 
direct estimate of the errors introduced on 
terminating it at any given term. 

With regard to the large values of p, the 
situation is still more favorable. For in the 
limiting cases © and o2/o0,=0.01 the use 
of the Hankel functions gives high accuracies 
with only two or three terms of the series. In 
the intermediate cases for o2/¢,;<1 the series of 
Eq. (15a) for the smaller values of p(p~2) begin 
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to diverge after the first two terms and must 
necessarily be dropped thereafter; since it may 
be proved that the remainders of the asymptotic 
expansions developed here are smaller in absolute 
value than the last term retained, the accuracy 
is not impaired (for o2/0,:=0.20, p=2, the 
remainder after the first two terms was 0.1 
percent). When o2/a,>1 the series do not diverge 
so rapidly and even though @p<1, three terms 
in Eq. (16) will in most cases already leave 
errors less than 3 percent. 

These results compare favorably with the 
laborious calculations required when the series 
of images are summed directly. For the con- 
vergence of these series, except when k~0, is so 
slow that to obtain accurate results as many as 
100 terms may at times be necessary.” 


II. PoTENTIAL DISTRIBUTION IN A STRATIFIED 
EARTH WITH MorRE THAN Two LAYERS 


A. General theory 


When there are more than two layers to be 
considered the numerical discussion immediately 
becomes much more complicated. The formal 
analytical expressions, however, may be readily 
extended from the two layer case treated above. 
To derive first these general formulae we may 
suppose that the potential distribution on the 
earth’s surface has again been found as an 
integral over Bessel functions and may be 
expressed in the units of Eq. (1a) as: 


b=1/e+2f fla) J(2pa)da, (19) 


where f(a), in the case of an n-layer earth is the 
ratio of two determinants of order +1 and 
contains as parameters the conductivities and 
depths of the m layers. For large values of a, 
f(a)~e-**, which insures the convergence of the 
integral in Eq. (19). 

We may now obtain the formal expansion of 
#—1/p as a power series in p by again replacing 
Ji(2pa) by its power series. We have then as 
the expansion to be used for small values of p: 


13 Cf. Trans. Amer. Inst. Min. Met. Engr., “Geophysical 
Prospecting,”’ 1932, p. 440. 
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1 (—1)"p" 
2mf(a)da. (2 


For large values of p we develop a general 
theorem for the asymptotic expansion of integrals 
of the type occurring in Eq. (19), which is a 
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where a=a, are simple poles of f(a) with 

residues b,. Eq. (19) then takes the form: 
Jo(2pa)d 


We now observe that: 


natural extension of the particular development 
used in the case of the two layer earth. Thus we 
suppose that f(a) can be expanded in a series of 
partial fractions in the form: 


f(a) =f(0) + {1/(a—an)+1/an}, 
+2 Io(—2pan) — Yo(—2pan)}. (23) 


Since we are interested mainly in the asymptotic expansion of ® with respect to p we may apply 
again Eq. (13) to obtain: 


Ji(2pa)da mr _ 


a+k 


where Hy is again the Struve function we met 


in Eq. (10a). becomes, then: 


1 1 
p 0 (— pa,)**! 


Now from Eq. (21), it follows that for m>0, f?"(0) =I'\(2m+1)>°b,/(—a,)?"*. 
We therefore have finally: 


+0(p- 


1 1» —1)"F?(m+4)f2"(0 
p 0 I'(2m+1) 


(24) 


Since Eq. (24) involves only differentiations of the function f it will not be difficult to apply even 
to complex forms for f. Eq. (20), on the other hand, involves the evaluation of a series of definite 
integrals which in general may be as difficult to evaluate as the single integral of Eq. (19). For 
small values of p it will therefore be necessary to generalize the method used in the case of the two 
layer earth. Although this method was essentially equivalent to summing the power series expansions 
of the individual image potentials, it is more convenient to derive the general expressions for the 
terms than to use the first few series of images already listed by Hummel,’ for a three layer earth. 


B. The three layer earth 
For a three layer earth the function f may be written in the form™ 


J (a) ; 
agit, A 


where: 


Ni =1—kee™; Ay — 1+ 


(25) 


h=he/hy>1.) 


The general evaluation of the coefficients of p’™ in Eq. (20) with the f(a) of Eq. (25) is 
unfortunately rather difficult. One can in a formal manner expand f(a) as: 


2(Not+k3N:1) Ai\" 
0 Ao 


4 Cf, Whittaker and Watson, Modern Analysis, p. 134. In the appendix we shall give a more general proof applying 
not only to functions f(a) as e~®* which do not permit the partial fraction expansion of Eq. (21), but also the 
expansions corresponding to Bessel functions of the second and third kinds and to functions of non-zero order. 

% Stefanesco and Schumberger, reference 4. 
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in which the term independent of k; will give 
the potential due to a two layer earth (k;=0, 
o3=02). The higher order terms can then be 
expressed as derivatives of this two layer 
potential with respect to kz. In this way one 
may obtain a power series expansion in p of the 
series of Hummel. The coefficients of this series, 
however, will themselves be double series with 
respect to powers of k, and k;, with unsatis- 
factory convergence characteristics. It will there- 
fore be more convenient to treat the potential 


for small values of p separately for special 
integral values of h, but by a method which not 
only leads to more satisfactory series but which 
also’ permits an extension to rational values 
of h. 

On the other hand, with regard to the series 
of Eq. (24) for large values of p, the general case 
offers no difficulties in principle, although the 
successive differentiations of f(a) become rather 
lengthy. Thus, putting in the value of f(a) as 
given by Eq. (25) we have for the first two terms: 


Pp=— 


93 


(2? — 03”) {01 — 02+ 02h}? — 02”) 
| (26) 
402°03"p? 


where it is to be recalled that h=he/h. 

As we should expect, if o2=03 or o1:=02, Eq. 
(26) reduces to Eq. (16) which we derived for 
the two layer earth. For very large distances 
®p approaches the value o;/0; and is independent 
of the conductivity of the second layer. Thus if 
o\= 03, Pp has, for large p, the same value as in 
a homogeneous earth. Its approach, however, to 
the value 1 does depend on a2, coming down 
from higher values if o2<o;, and up from lower 
values if o2>0;. On the other hand, for smaller 
distances the second layer does come into play 


and its effect increases rapidly with its thickness, 
as is clear from the second term in Eq. (26). 
We shall now return to the potential distribu- 


tion for small values of p for two special values 
of h. 


Power series expansions for special values of h. 


h=2. When the thickness of the second layer 
is exactly equal to that of the upper layer the f(a) 
as given by Eq. (25) can be re-expressed in a 
form which permits an application of the method 
used for a two layer earth. Thus we have 


f(a) = 


k3e~** +- koe* 


where: 


= — ke(1 k3) + [ke?(1 k3)?+4k; }} = 2k3/s1, 
2k3e°2 = — k3) [ke?(1 k3)?+4k; = 2k3/Se. 


e~ v2 


(27) 
sinh (a+71) sinh 


(27a) 


Carrying through now the same expansions and integrations as in Section A of Part I, it is found, 


by means of Eq. (19), that for the present case: 


1 
b=-+ 


ke k; 


+ E log 
p [ko?(1—k3)?+4k; }} 1—s; Ay 


2 ks ks 
+— log (1—s2) -— log (1—s)) 
(28) 


4 


If ks=0(c3=02), Eq. (28) easily reduces to the ® of Eq. (5) for a two layer earth; and if 
k2=0(¢2= 01), Eq. (5) is again found with p changed to p/2 and a factor of 2 in ®. 

If the lower layer has an infinite resistivity (0;=0; ks= +1), Eq. (27) gives an infinite constant 
in addition to the series, as is to be expected since f(a) has a pole at a=0 if kz=1. Eq. (28) reduces 


| 
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then to: 


1 
(k;=1, o3;=0) =-+ log 2—(ke+1) lim log (1 — 51) 
p 


dn ( 


In using these formulae the procedure followed 
in Part I may be applied in separating out and 
summing the terms for »=1, and 2, thus in- 
creasing the radius of convergence of the series. 
Since |s;| and |s2| are less than 1, the same 
function f,,(k) occurring in Eq. (6) can be used 
with Eq. (28) after separating out the n=1, 2 
terms, except that the arguments will be s; and 
S2 instead of the k of Eq. (6). 
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28. 
(28a) 


1 


For large values of p, the more general Eq. 
(26) may be used with h set equal to 2. 

h=3. When the second layer is two times as 
thick as the upper layer the analysis can again 
be carried through in detail, although it will be 
convenient to make it somewhat more formal 
than in the case of h=2. Thus we shall suppose 
that the denominator of f(a) has been factored, 
so that f(a) can be rewritten as: 


— 


f(a) ~ (e~% e771) e772) e773) 
(se— 51) (S3—51) lt (S1—S2) (Ss—S2) $3) (S2— 
sinh (a+71) sinh (a+72) sinh (a+72) 
where 


e 


= $9; 


= S3. 


Again applying the method of Section A, Part I we find: 


1 2 
p 0 


—1)"cnp™ 


ke ks 
(n+ 


(30) 


— Se) — 53) }. 


Here, too, the constant terms for m=0 and 
those for n= 1, 2, may be taken out and summed 


electrode 


| 
h, a) Oo, 
h, 
(2) 
0, 
v. xy 
Z(w= %n) 


Fic. 4. Diagrammatic representation of a three layer earth. 


separately and the remaining sums with respect 
to expressed in terms of f,,(s). Although 52, 53 
cannot be expressed in simple explicit forms, 
they may be readily obtained in any particular 
case as the reciprocals of the roots of the tri- 
nomial denominator of f(a) set equal to zero. 
And again ® for large p may be found by means 
of Eq. (26) with h set equal to 3. 

It may be noted that insofar as the roots of 
any polynomial may be found graphically the 
above method may be extended directly to 
systems where / is any rational fraction, although 
it will of course become increasingly complex as 
the degree of the polynomial increases. 


Lower layer of zero conductivity 


_ When the lower is of zero conductivity (o;=0, 
k3= +1) the potentials diverge, as in the case 


of the two layer earth. This may be seen on the 
one hand from the asymptotic expansion of Eq. 
(26), which becomes infinite when o3;=0, or from 
the form of f(a) of Eq. (25), which has a pole at 
a=0 for +1. 

For small values of p we have seen how the 
divergence of ®@ may be separated out into a 
logarithmic term. Thus for h= 2, Eq. (28a) gives 
the infinite part of as —(1+2) lim log (1—s}). 


And similarly in the case h= 3, it may be shown 
by summing in Eq. (30) the terms with m=0 
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—2(1+ke) 


lim 1 1—s)). 
ae" og ( 1 


For large values of » we must find a new 
asymptotic expansion, as Eq. (26) diverges 
everywhere for o;=0. The situation is similar 
to that encountered in the case of Eq. (16) for 
a two layer earth, and may be treated in the 
same way. 

Thus, for h=2 we see from Eq. (25) that the 
poles of the two parts of f(a) are: 


and taking s, as the reciprocal of the root of An (31) 
value unity, that the infinite term of is: with residues: 
(—1)"e-" — (— 1) + 
; (31a) 


}! 


[ko?(1—ks)?+4ks }! 


Noting that for ks=1, y:1=0, y2= 77/2, and that for each set, b, = b_,= do, it is readily found that: 


lim 
(k3—>1) 


{1+f(0) + =0. 


Applying now Eqs. (23) and (17) we have finally, 


lim, lim) log 2y,;—0.5772+log 2/0} | 


+23. nKo(nxp) —1)"Ko(n7p). 


(32) 


For h= 3, we return to Eq. (29) and note that for k3=1: 


1=0; y2=—10/2; 


ba(vi) = (1 /(3— ke); 
so that, again: 


¥3=10/2; 
ba(y2) = = (1 — ke) /(3 — ke) 


cos 6 = (ke—1)/2 (33) 


fim, =0. 


Applying once more Eqs. (23) and (17) we finally find that 


o3->0 


| + o(2nxp—Op)} 


lim (h=3) =[2(1+h2)/(3— ke) ] {—lim log 2y,—0.5772+log 
1 


(34) 


Similarly, for h=4, it may be shown without difficulty that for o,=0(k;=1) the poles of f(a) 


and their residues are: 
a,=n7ri+0; 


—7i/2; 


— 10/2; 
1—k,? ko}. 


nxi+i0/2 
(35) 


cos 
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Putting these in Eq. (23) and making use of Eq. (17) we find that: 


jim) =4) ] —lim, log 2y,—0.5772+log 2/p 


Before presenting the analogous formulae for 
the case of a lower layer of infinite conductivity 
(o;= ©, k;= —1) it may be observed that Eqs. 


convergent and involve no logarithmic terms or 
infinite constants. The method of derivation, 
however, is identical with that for o,=0 and 


(32), (34) and (36) all reduce, as they should, will be only briefly indicated. 
to Eq. (18), derived for a two layer earth, when Thus for h=2, the poles and residues of f(a) 


ke is set equal to 1(¢2=0), or O(o2=01). | are: 
Lower layer of infinite conductivity |. (37) 
bn=3; 


When the lowest layer is of zero resistance 
closed series of Hankel functions may again be By using once more the general series of Eq. 
obtained. As in the corresponding case fora two (23) and reducing by means of Eq. (17), it is 
layer earth (cf. Eq. (11c)), the series are strictly found that: 


0; h=2) (38) 
0 1 


For h=3, we have similarly: 


Qn =(2n+1)xi/2; nri—i0/2; nri+i0/2; (0) 


so that: 


4(1—ke), 4(1+he) 
~, h=3)= o((2n+1)rp)+ nK (40) 
3+hk. ° | 


Finally, for the case h=4, we have: 


cos = {ke+ (ko? +8)'}/4; cos ke— (ke? +8)!} /4, 
_ (he? +8)! _ Shot (he? +8)! (he? +8)! (he? +8)! 
+8)! 4(ko? +8)! 4(ko? +8)! 
so that: 1+/(0)+),/a,=0 


v 


(41) 


ar 


| 
‘ a 
0 
li 
g 
a 
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and 
= co: => = _ 
> (42) 


— 3ko+(k2?+8)! | 
+ > nKo(2nrp+6ep) + nK o(2nrp—Oep) +. 


(ke? +8)! 


1 


Again it may be observed that Eqs. (38), (40) and (42) reduce to that for the two layer potential 
(cf. e.g., (11c)) both when ky is set equal to —1(¢2= ©) or O(¢2=01). 


Numerical treatment 


To illustrate the analytical results obtained 
above the computations for a number of cases for 
which numerical results have as yet not been 
given in the literature have been carried through 
in detail and will be presented in graphical form. 

In Fig. 5 curves are given for &p for the cases 
of a slab of thickness equal to that of the top 
layer (k=2) immersed in an otherwise homo- 
geneous earth (¢;=o;). The two cases are for 
and o2=5e0;. As already suggested in 
our discussion of Eq. (26), ®p for these cases 
approaches 1 for large p, with #p>1 for 


and @p<1 for o2.=5e;. For small p 
Eq. (28) was used with the terms for n=1, 2 
summed separately. In the case ¢2=50,, s; and 
Sg were real and the resulting sums for m were 
expressed as f,,(s). In the case ;= 51, Se were 
complex but the series reduced quite simply, 
although the f,,(s) could not be used. In both 
cases only four terms in the series with respect 
to m gave high accuracies even for p=2.5. For 
large p the extensions of Eq. (26) were used. 
The higher derivatives indicated in Eq. (24) 
were most easily derived by expanding separately 
the numerators and denominators of f(a) in 


oat \ 


Fic. 5. Surface potentials for the case of a slab of thickness equal to that of the top layer, immersed in 
an otherwise homogeneous earth. 


144 


MORRIS MUSKAT 


>= 
oN 
P ~ 
7 
as 10 20 25 30 40 50 


Fic. 6. Surface potentials for a three la 
and the conductivity ratios of the two up 
h:—h, =thickness of second layer; p = (radial distance from electrode) /2/. 
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Fic. 7. Surface potentials for a three layer earth in which the lowest layer has a zero conductivity and 
the conductivity ratios of the two upper layers are 1: 5 and 5: 1. 
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powers of a and then getting the coefficients of 
2” in the ratio by well-known algebraic rules. 
In the case of o2.=50,, the convergence of the 
asymptotic expansion was made very rapid by 
using the analogue of Eq. (15a), i.e.: 


1 
Lo(—2pao) Yo(—2pao)} 
1» 


43 
=m I'(2m+1) p?"*! 


bol (2m+1) 


a 


where the notation is that of Eq. (23), ao being 
the pole of f(a) of smallest absolute value. 

Cases in which the lower layer is of zero con- 
ductivity (¢;=0) are illustrated in Figs. 6 and 7, 
for middle layers of thickness, 1, 2 and 3 times that 
of the upper layer. Eqs. (32)—(36) were used in 
computing these. The inflections in the curves 
for ¢2>0, and the almost exact overlapping of 
all the curves for p<1 may be noted. Of more 
interest, however, from the analytical point of 
view is the fact that the equations (32)—(36) 
(as well as Eqs. (38)—(42) for o;= ©) are rigorous 
expressions for ® and furthermore converge so 
rapidly that even for p=0.25 less than eight 
terms in the series over Ko will give errors less 
than 0.05 percent. For p=1 only two terms give 
equally high accuracies. As, by definition, p= 1 
for p=0, these series were very satisfactory for 
numerical treatment over the complete range 
of p. 

Similar curves for o;= © are given in Figs. 8 
and 9. Of particular interest here is the fact that 
although all the curves approach 0 as p becomes 
large the curves for o2>0; lie much closer 
together, just as in the case of the two layer 
earth (cf. Fig. 3). As to the convergence of the 
series, they are just as satisfactory as those for 
o3;=0. 

As a single generalization of the detailed 
analysis developed here for two and three layer 
earths, we shall simply mention that since for 
an layered earth the po- 
tential, by Eq. (24), for an m layered earth will 


approach the value 
&=0,/onp (44) 


as p becomes very large. Thus, in general, as 
already suggested by others, the potential at 


145 


large distances from the electrode is determined 
only by the topmost and deepest layers. And 
finally it may be noted that the asymptotic 
expansions developed here can be applied to 
continuously varying conductivities as well as 
those varying discontinuously. 

In concluding, the author wishes to mention 
that after all the analysis of the present paper 
had been developed and written up, and while 
the author was waiting for the completion of the 
numerical work, a paper by L. V. King on the 
same problem appeared in the February Ist 
number of the Royal Society Proceedings. 
Although King’s analysis in Sections 7-11 is 
almost identical in many respects with that 
developed by the author there are after all 
enough differences in the final results that it 
has seemed best to leave the author’s original 
manuscript unchanged. Among the results de- 
veloped by the author and not given by King 
are: (1), detailed discussion of the power series 
for small values of p (Eqs. (6), (20), (28), (30)); 
(2), expressions valid when the electrode is not 
on the surface and for potentials below the 
surface (Eqs. (8), (11), (12)); (3), the closed 
expansion for ® in terms of Struve’s function, 
which is very useful in getting the potentials for 
o3=0, ©, as well as for evaluating, in general, 
a large class of Bessel function integrals (Eq. 
(23)); (4), explicit expansions for the potential 
for a three layer earth (Eqs. (34)-(43)); (5), 
explicit expressions for the asymptotic expan- 
sions in terms of the conductivities (Eqs. (16), 
(26)); (6), thirty numerical cases represented 
in graphical form; and (7), a general proof of 
the expansion of Eq. (24) which permits the 
method to be applied to continuously varying 
conductivities (Appendix). On the other hand 
it should at the same time be mentioned that 
King does include discussions of field graphs and 
alternating current problems, which have not 
been considered in the present paper. 

The writer is indebted to Professor Harry 
Bateman for helpful suggestions with respect to 
some of his earlier work on this problem, to 
Mr. R. D. Wyckoff, Dr. P. D. Foote and Dr. 
E. A. Eckhardt for permission to publish this 
work, and to Mr. M. W. Meres for the pains- 
taking care with which he has computed and 
plotted the data for the various figures of this 
paper. 
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Fic. 8. Surface potentials for a three layer earth in which the lowest layer has an infinite conductivity 
and the conductivity ratios of the two upper layers are 1 : 2 and 2: 1. 
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Fic. 9. Surface potentials for a three layer earth in which the lowest layer has an infinite conductivity 
and the conductivity ratios of the two upper layers are 1: 5 and 5: 1. 
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APPENDIX 


As may be readily verified by trial, Eq. (24) is valid 
not only for such functions f(a) as permit an expansion in 
partial fractions, but also for such functions as e~*, e~6* 


or Jy(a). That it has a more general validity—although 
Eq. (23) applies only to functions with partial fraction 
expansions—may be proved as follows. 


By Watson, p. 170, we have for: —}<R(v)<} 


Pip) = (ap) f(a)da= sin aplda 
Now the second integral has the asymptotic expansion and 
(successive integrations by parts): =. f"(0) 
J, sin aptda= (46) we get: 
Putting this in P,(p) and integrating with respect to 
we find: 
» @ rip’ & 1) /2—v)fr(0) 


which for »=0 is identical with Eq. (24), taking into 
account the difference between the definitions of p and 
f(a), here and in Eq. (24). 

Noting similarly that for —}<R(v)<} 


cos apt 


a” Y,(ap) (2 1)** j 


and that: 


0 


(- 


we get: 


= 1) 


2” 
wip i) 
And finally, since for —}<R(v) <1 


T(3)p” ~apt(y2— 


(48) 


When f(a) is odd or even Eqs. (47) and (48), respectively, 
break down. As already seen in the text this is due to 
the fact that in these cases the integrals P and Q have 
exponentially decreasing asymptotic expansions. This may 
again be proved generally by considering the functions 


sin apt |. 


as f(a) is odd or even. Thus, if f(a) is odd and f(a)~ 1/a"*® 
for |a|~ x acontour integration of the integral shows that 


where a, are the poles of f(a) in the first quadrant, 8, 
the poles in the 4th quadrant and e,=2 if the real parts 
of (@,, 8.) are positive and e,=1 if they are 0. It follows 
at once that J vanishes exponentially with large p, and 
similarly for the integral with cos apt if f(a) is even. 
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The Crystal Structure of Uranium 


Tuomas A. WiLson, Union College, Schenectady, N. Y. 
(Received February 6, 1933) 


A method capable of solving x-ray powder photographs 
derived from crystal structures of low symmetry has been 
described. By the use of this method, the crystal structure 
of uranium has been obtained. Uranium crystallizes 
according to space group C2,* of the monoclinic system, 
with 2 atoms in each crystallographic cell. The dimensions 


of the cell are: a=2.829A, b=4.887A, c=3.308A. Axial 
ratios are: 0.5791 :1:0.6771. The x-ray density of 
uranium is 19.05. The uranium atoms appear to be prolate 
spheroids with 3.425A as the major axis and 2.826A as 
the minor axis. 


HE crystal structure of uranium has been 

reported by McLennan and McKay! as 
body-centered cubic, with a=3.43 Angstrom 
units. Knaggs, Karlik and Elam? report that a 
better agreement with the experimental data is 
afforded with a=5.50A. In Wyckoff’s compila- 
tion’ of crystal structure data uranium is noted 
as being not cubic. The spacings observed by 
McLennan and McKay are in very poor agree- 
ment with either of the cubic structures sug- 
gested, and their reported spacings do not agree 
with those observed on his x-ray films by the 
author. Their procedure of preparing the speci- 
men for analysis is open to question. 

Uranium is representative of a large class of 
substances possessing lower crystal symmetry 
of which single crystals are either difficult or 
impossible of preparation. Hence the crystal 
structure must be determined from an analysis 
of x-ray photographs obtained by some variety 
of the powder method. Runge‘ has proposed 
such a method, which seems to have certain 
objections preventing its wide adoption. Another 
method of performing this analysis, and the 
results obtained by it, form the subjects of this 
paper. 


METHOD OF ANALYSIS OF THE 
POWDER PHOTOGRAPHS 


The method involves the resolution of the 
simple atomic planes of a crystal structure into 


1 McLennan and McKay, Trans. Roy. Soc. (Canada) 
C111] 24, 1 (1931). 

? Knaggs, Karlick and Elam, Tables of Cubic Crystal 
Structure, London, 1932. 


the [100], [010] and [001] zonal families, and 
determining the interplaner relationships neces- 
sary in each family to satisfy the observed x-ray 
reflections. 

Consider the nearest equivalent points of some 
crystal structure, with origin at 0. These points 
define a parallelopiped called the ‘‘simple cell” 
in this paper. In a homogeneous space lattice, 
such as the uranium lattice, such a cell contains 
one equivalent point with the coordinates 000. 
Hence the structure vector becomes identical 
with the “loading factor’ of Niggli,> and no 
extinctions are produced. 

Let D,, Do, --- D, be the interplaner spacings 
obtained for this simple cell by x-ray powder 
diffraction methods. Then 


D? : D,? : D3: +--+ D,? 


where F,, F2, F;, --- F, represent the F vectors 
of the reciprocal lattice. By definition, 


F= K,,h?+ Kook? + hk 
+2Ki3hl+2Ko3kl (1) 


where h, k and / denote the plane indices. The 
smallest values of F will generally be obtained 
when F;, F, and F;=Ki:, Kez and K33, respec- 
tively. Hence the largest values of D on the x-ray 
film will correspond to the face plane spacings 
of the cell. 


3 Wyckoff, Int. Crit. Tab. I, 341 (1926). 

* Runge, Phys. Zeits. 18, 509 (1917). 

5 Niggli, Handbuch der experimental Physik VII, Part 1, 
239 (1928). 
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CRYSTAL STRUCTURE OF URANIUM 


Let e= D;/D,=1/F;! =f= D2/D,=1/F,' =g 
= D,/D,=1=1/F;' define the three vectors e, f 
and g from 0 which represent the face plane 
perpendiculars of a cell of such size that the 
longest perpendicular distance is regarded as 
unity. (Fig. 1a.) These vectors define a reduced 
perpendicular lattice.6 In such a lattice, let 
d[100], d[010] and d[001] represent the inter- 
planer spacings of the [100], [010] and [001 ] 
zones, respectively. These quantities will be 
defined by the equations: 


1/d?[100]= +P+ 2Ko;3'kl (2) 
1 /d?(010 J= +P +2K,;/hl (3) 
1 /d?(001 (4) 


in which the terms 2Ko;’, 2K,3;’ and 2K,’ are 
still to be determined. Let 6 and e refer to the 
angles shown in the ef plane of Fig. 1a. Then 


2K = sin? 5/sin® e— (5) 
and if 1/m?, 
tan e=m sin 6/(1—m cos 4). (6) 


Similar relations exist between 2K,;’, Ki’, B 
and 6; and 2Ko;’, Koo’, y and a. 


O<e 
We 


b 
Fic. 1. Reduced lattice of uranium. 


Having obtained Ko’ in Eq. (2) in the manner 
described, the value of 2K; can be obtained if 
some observed x-ray reflection on the powder 
photograph is considered as showing the spacing 
of the (011) plane, equivalent to d(011) of the 
reduced lattice. This reflection may occasionally 
be one previously used in connection with the 
face planes. A short series of increasing F- 
functions is calculated, and from these the 


® Reference 5, p. 108. 
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corresponding series of decreasing D-functions. 
If this series agrees with observed interplaner 
spacings on the x-ray film, the series is correct. 
If such agreement is not obtained, another 
observed x-ray reflection is considered as D(011) 
and another equation and series are set up. 
Usually only a trial or two is necessary. 
Identical processes are carried out for Eqs. 
(3) and (4). When these equations have been 


found, the general equation of the simple cell 
becomes 


(7) 


As a final check, the interplaner spacings of the 
other zone families are calculated and compared 
with the observed x-ray spacings. 

If V represents the volume of the simple cell, 
the expression 


V=D{e-(fXg)] (8) 


gives the volume allotted to each atom in the 
homogeneous space lattice. 


PREPARATION OF THE SPECIMEN 


The uranium used in this work was in the 
form of a button. It had been prepared by the 
method of Goggin, Cronin, Fogg and James.’ 
Since filings of uranium catch fire in the air, the 
x-ray photographs were obtained by reflection 
from a plane surface of the metal, formed by 
filing the uranium flat in a hydrogen atmosphere. 
A bright silvery surface was thus obtained. This 
surface was not coated to protect it from the 
oxidizing action of the air. The exposures were 
obtained in 23 hours. No visible tarnishing of 
the metal surface occurred during this period. 
For additional checks, to show absence of 
sufficient oxidizing action to affect the x-ray 
data, photographs were obtained from powder 
specimens of the uranium oxides supplied by 
different chemical firms. Films were also prepared 
from filings of the metal which had been allowed 
to oxidize in the air at ordinary room tempera- 
ture. Other films were made from filings heated 
in the air on a hot plate. None of these films 
gave x-ray patterns which agreed in any way 


7 Goggin, Cronin, Fogg and James, Ind. Eng. Chem. 
18, 114 (1926). 
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with those obtained from the filed surface of the 
metal itself. 

The photographs were taken by Dr. W. P. 
Jesse of the General Electric Research Labora- 
tory. Unfiltered iron radiation was used from a 
tube operated at 16 m.a. and 38 kv. The radius 
of the Pye camera was 45.7 mm. The filed 
surface of the block of uranium was centered in 


the camera as precisely as possible and was 
rotated by hand during the exposures. 


X-Ray DaTA 


Table I contains the readings made from the 
second of two exposures. The lines were some- 
what clearer on this film, although the lines 
recorded on the first film agreed in position with 


TABLE I, X-ray data. Derivation of the simple uranium cell. 


(1) (2) (3) (4) (5) (6) 


Readings of film, mm 


(7) (8) (9) (10) 


Observed Observed Cale. 


from zero beam Intensity Radia- D-values  d-values d-values 
Left Right Av. estimated __ tion A A Equation Atomic planes 
oad 32.0 32.2 w B 2.538 
33.8 33.3 33.5 w B 2.445 
35.9 35.8 35.8 ms a 2.535 1.0000 1.0000 10, 11 (100), (010), (001) 
36.3 36.3 36.3 w B 2.261 
37.2 37.1 37.2 ms a 2.443 0.9637 0.9637 10 (110) 
40.6 40.4 40.5 s a 2.252 » 0.8884 0.8884 11 (101) 
47.3 47.3 47.3 w B 1.774 
52.9 53.0 53.0 s a 1.764 0.6959 
56.0 —- 56.0 f B 1.524 
60.5 — 605 vf 1.427 
63.0 62.9 63.0 nis a 1.519 0.5992 0.6019 sogye 11 (701), (10T) 
1.378 0.5949 10 (110); (110) 
68.0 68.0 68.0 w a 1.427 0.5629 0.5631 10 (120) 
70.9 —- 70.9 w a 1.381 0.5448 0.5361 11 102) 
82.0 82.0 82.0 ms a 1.236 0.4876 0.5000 ho 4909* 10, 11 (200), (020) 
0.4818 f 2 (220), (002) 
91.1 — 91.1 w a 1.151 0.4540 0.4490 11 (202) 
96.7 a 1.109 0.4375 
B 1.006 
117.0 117.0 f a 1.008 0.3976 0.3926 11 (201) 


s—strong, ms—medium strong, w—weak, f—faint, vi—very faint. Average values shown by *. 


those on the second. Readings were made of the 
positions both right and left of the zero beam 
as shown in columns 1 and 2 of the table, and 
the average of the readings was taken as the 
basis of the calculations (column 3). The wave- 
lengths of the Fe Ka, and the Fe K@& radiations 
(1.9323 and 1.7527A, respectively) are those 
given in Compton.’ The D and the d values 
calculated from the x-ray data are tabulated in 
columns 6 and 7. 

The small number of interplaner spacings 
shown on the film indicate a rather high sym- 
metry in the simple uranium cell. Such symmetry 
would be imparted to the ef plane (Fig. 1a) 
if e=f, or if 6=«. From previous work on the 
series of x-ray reflections observed under each 
of these conditions it seemed more probable 
that a series in which e=f would fit the observed 


Compton, X-rays and Electrons, p. 396, New York, 
1928. 


data. The general equation for the interplaner 
spacings of the [001] zone then becomes: 


1/d?(001 J=h? +k? + 2K (9) 


Considering d(110)=0.9637 in the reduced lat- 
tice, 2Ki2’= —0.9232, hence Eq. (9) becomes 


(10) 


Comparison of the interplaner spacings calcu- 
lated by means of Eq. (10) and those observed 
on the x-ray film is shown in columns 7 and 8 
Table I. The atomic planes of the simple cell 
are listed in column 10. Since in this equation 
K,)'= 2K 12’=sin? 6/sin? e—2, 2e= 180 
— 6, from which 6= 62° 31’. 

By applying the same reasoning to the [010] 
zone, with the value of d(010)=0.8884, 2K; 
— 0.7576, and Eq. (3) becomes 


1 /d?(010]= — 0.757601. (11) 
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TABLE II. Comparison of the observed and calculated D values of uranium. 


(1) _ * (2) (3) (4) (5) (6) 
Calculated values of simple cell Observed values Calculated values of monoclinic cell 
Interplaner Interplaner 
spacings (D), Spacings Inten- spacings (D) 
Planes in A in sity Planes in A 
(010) (100), 2.535 2.535 ms (011), \ 2.535 
(110) 2.443 2.443 w (020) 2.443 
(101), (011) 2.252 2.252 s (111), (171) 2.281 
(112) 1.772 1.764 s (120), (120), 1.760 
(221) 1.718 20) . 
(101), (O11) 1.526 1.519 ms iy, (171) 1.538 1.510* 
(210), (120) 1.427 1.427 w (30, (031) 1.435 ie 
(102) 1.433} 
(102), (201), } 1.381 1.381 w (131), (131), 1.383 
(021), (012) ’ (211), (122), } 1.359 | 1-368* 
(020), (200) } 1.289 
’ , 1 . 
(002) 200), (022 1.270* 
(2) 1.236 ms } 1.268 
(202), (022) 1.137 1.151 w (130) (131) 1147 \ 
(123) 1.109 (220), } 1.125 
ab. (211 1.000 102 1.086 
(102), (301) 0,993 9-996" 1.008 f 113%, (113) | 1.08% 
(012), (021) (231), (231) 1.071 


s—strong, ms—medium strong, w—weak, f—faint. Weighted mean values indicated by *. 


Results obtained from this equation are likewise 
tabulated in column 8. The calculated value of 


B is 67° 43’, 


2.226 AU 


Fic. 2. The uranium cell. 


So many of the reflections are now accounted 
for that it seems probable that 2K,;=2Ko;, 
hence B=. So the projection of g on the ef 
plane coincides with the bisector of 6. This 
situation produces the symmetry expected in 
the simple uranium cell, although the cell is 
triclinic. 

Since D,=2.535A from the x-ray data, a 
rearranged form of Eq. (7) for the simple 
triclinic uranium cell may be written: 


2.535?/D?= h? +k? +P —0.9232hk 
—0.7576(h+k)l. (12) 


Columns 2 and 3 of Table IT show the agreement 
obtained between the observed D values and 
those calculated by Eq. (12) for the planes 
shown in column 1. 

The cell thus derived is shown in Fig. 1b. 
The volume of this simple cell (Eq. (8)) is 
20.53A%. The calculated density of uranium is 
19.05, compared with the usually reported 
density of 18.68; a difference of 2.2 percent. 


THE CRYSTALLOGRAPHIC CELL OF URANIUM 


By a change in the choice of coordinate axes, 


| 
= 


152 THOMAS A. WILSON 


this simple cell may be converted into the 
monoclinic cell shown in Fig. 2. The crystallo- 
graphic axes, a, b and c, are indicated. The 
face plane perpendiculars of this cell are EK, 
Ell and EL, with the lengths 2.535, 2.966 and 
4.887A, respectively. 5= 63° 26’. The equation of 
this cell is: 


2.966? /D® = 1.368h? +0.368k? +1? —1.046h1. (13) 


The coordinates of the equivalent points of the 
cell are 000, 033, which locates uranium in space 
group C2 of Wyckoff's tabulation.’ According 
to the structure vector for this configuration, 
the planes producing reflection are those in 
which h, k and / satisfy the following conditions: 


hki all even, or all odd; 
h even, ki odd; h odd, ki even. 


Columns 3 and 6 of Table II show how the 
interplaner spacings of the permitted planes 
listed in column 5 agree with the observed 
spacings. The number of cooperating planes 
producing reflection is in satisfactory agreement 
with the observed intensities, if allowance js 
made for the natural decrease of x-ray intensity 
with increasing angle of reflection. 

In this cell, a= 2.829A; b= 4.887A; c= 3.308A, 
The axial ratios are therefore 0.5791 : 1 : 0.6771. 

As shown in Fig. 2, this structure corresponds 
to a close packing of prolate spheroids with 
minor axes= a= 2.826A, and major axis= 2E///3! 
= 3.425A. This appears to be the probable shape 
of the uranium atoms in the lattice. 


* Wyckoff, The Analytical Expression of the Results of 
the Theory of Space Groups, Washington, 1922. 
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Magnetic Properties of the Iron-Nickel Alloys Under Hydrostatic Pressure 
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A series of well annealed iron-nickel alloys in nominal 
steps of 10 percent Ni, including pure iron and pure 
nickel, has been examined for the effect of pressures, up 
to 12,500 kg/cm’, on the magnetic flux density. The 
application of pressure to a well annealed specimen causes 
a relatively large change in flux, usually a decrease and 
an incomplete recovery on release of pressure. The recovery 
is reversible in the sense that the same curve of flux change 


is traced with each series of pressure applications, provided 
the measuring pressure does not exceed that first applied. 
The 30 percent Ni alloy is rendered nearly nonmagnetic 
by a pressure of 12,000 kg/cm*. In the remaining metals 
pure iron exhibits the greatest and 90 percent Ni the 
smallest pressure change of flux. At constant field the 
different alloys exhibit various types of pressure changes 
of flux; linear, nonlinear and types exhibiting hysteresis. 


METHOD AND EXPERIMENTAL PROCEDURE 


HE high-pressure technique of Bridgman! 

has already been applied to the examina- 

tion of the ferromagnetic properties of the pure 
metals, iron, nickel and cobalt, by Dr. Chi-Sun 
Yeh.2 The method of the present investigation 
differs from that of Yeh principally in an 
improved annealing treatment. Small rings of 
rectangular section with an average radial thick- 
ness 0.2 cm, axial width 0.5 cm and average 
outside diameter 2.0 cm were cut from a series 
of vacuum fused iron nickel alloy plates very 
kindly prepared for this investigation by Dr. 
L. W. McKeehan at the Bell Telephone Labora- 
tories. The alloys were prepared in nominal 
steps of 10 percent nickel from pure iron to pure 
nickel. The chemical analyses furnished by Dr. 
McKeehan with the alloys are given in Table I. 
After machining, each ring was heated to 
1100°C in a vacuum furnace, then slowly cooled 
to room temperature at the constant rate of 
30°C per hour. Secondary windings of 500 to 
800 turns of No. 40 silk covered wire were wound 
on the rings over an insulating layer of thin 
paper. The primary winding of 200 to 400 turns 
of No. 32 silk covered copper wire was placed 
on top of the secondary. The completed speci- 


1P. W. Bridgman, The Technique of High Pressure 
Experimenting, Amer. Acad. Arts and Sci. Proc. 49, 
627-643 (1914). 

2Chi-Sun Yeh, The Effect of Hydrostatic Pressure on the 
Magnetic Permeability of Iron, Cobalt and Nickel, Am. 
Acad. Arts and Sci. Proc. 60, 503-533 (1925). 


mens were then placed, one at a time, in a steel 
cylinder in which a pressure of 12,500 kg/cm? 
could be developed. 

The electrical measurements involve the fa- 
miliar ring ballistic method and need not be 
described here. For some specimens the compen- 
sating method of Yeh? was used, wherein the 
flux linkage between the two windings on the 
specimen was balanced by an equal and opposite 
linkage between two external air core coils. 


TABLE I. Chemical analyses of iron-nickel alloys. 


Specimen Nickel Iron Manganese 

numbers (%) (%) (%) 
M.N. 1, M.A.S. 0.07 _ 0.03 
M.A.S.R., M.A, 2 
M.A. 27 10.09 89.0 0.32 
M.A. 28 20.15 78.69 0.36 
M.A. 29 30.02 68.85 0.38 
M.A. 30 40.51 58.37 0.37 
M.A. 31 50.04 49.08 0.34 
M.A. 32 59.57 39.87 0.37 
M.A. 33 69.27 30.03 0.42 
M.A. 34 78.98 19.94 0.60 
M.A. 35 89.22 10.25 0.40 
M.A. 36 _ 0.30 0.42 


Two types of run were made, constant pressure 
and constant field. The first gave a B-II curve at 
a fixed pressure. In the second type the field 
was held fixed and the resulting change of B 
with pressure was measured. In general three 
constant pressure runs were made on each 
specimen at 30°C. A few were made at room 
temperature and a few at 75°C. The first run 
on each specimen, made at one atmosphere 
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immediately after annealing and called the 
primitive, serves to describe the original magnetic 
condition of the metal. The next run was made 
nominally at 12,000 kg/cm’. The third was 
made at one atmosphere after release of pressure, 
and served to discover the permanent effect of 
pressure which is quite large in nearly all the 
specimens examined. All ballistic measurements 
were made at constant reversal (i.e., the meas- 
uring field was reversed a number of times until 
the ballistic deflection became constant before 
taking a reading) except in a few cases where 
the first reversal was recorded. 

The effect of temperature and pressure in 
altering the flux and field measurements through 
changes in dimensions is small, not greater than 
0.5 percent, but the effect was corrected for. 


DISCUSSION OF RESULTS 


The tables of ballistic data and the detailed 
mechanical constants are too voluminous to 


0 /0 20 30 40 


include in this report. A thesis* presented at 
Harvard gives a full account of the work. The 
essential behavior of the alloys is clearly shown 
in the curve charts. The series of curves is not 
entirely complete. The 10 percent Ni alloy was 
omitted and the observations on the 70 percent 
Ni alloy were not completed because of lack of 
time. However, the essential behavior seems to 
be revealed clearly enough in spite of these 
omissions. 


Constant pressure runs 


In Figs. 1 to 8 are given the curves for the 
constant pressure runs. The general effect of a 
pressure of 12,000 kg/cm? is to cause a large 
initial decrease of flux density with only partial 
recovery on release of pressure. Pure iron, 60 
percent and 90 percent nickel alloys, and pure 


3R. L. Steinberger, Magnetic Properties of the Iron 
Nickel Alloys under Hydrostatic Pressure, Harvard Physics 
Thesis, February, 1932. 
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nickel, exhibit reversals in sign of the pressure 
effect—pure iron near JJ=10, 60 percent Ni 
between J7= 20 and J7=40, 90 percent Ni near 
H=1.5 and pure nickel near H=2. In the 30 
percent alloy the recovery is practically complete 
on release of pressure. This specimen is also 
unique because of the exceptionally large value 


at | 
lhe 
not 
ent 
of 
te 10000 
4 8000 
| 
4000 fj 
/ 
the 
fa 4000\— 
rge 
tial 
60 
ure H 
[ron 
“| 
200 
000 
6400 
iS 
| 
H 
0 
2 3. 
Fic. 4. 


156 


of the reversible change of flux with pressure 
which will be enlarged upon in the constant 
field run discussion. 
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Temperature effects 


The general effect of a rise in temperature in 
the range 30° to 75° is to cause a decrease in 
flux density. The droop in the constant pressure 
curves (Figs. 3 and 4) is due to the heating of the 
specimen by the magnetizing current in the 
restricted space of the pressure apparatus. In 
the 30 percent alloy the temperature coefficient 
is quite large, amounting to a decrease of 51 
maxwells/cm? per degree at /7=7.5, 30°C and 
at a pressure of one atmosphere. 


Constant field runs 


Figs. 9 to 18 list the constant field runs which 
reveal the percentage change in flux density at 
a fixed value of JJ as the pressure is raised and 
lowered. The constant field is not the same for 
all specimens. It was fixed at a figure near 
where the percentage change of flux is a maxi- 
mum (see Fig. 20). The constant field curves 
exhibit a wide variety of forms. They indicate 
both positive and negative changes of flux with 
pressure. Most of the curves exhibit a decrease 
of flux. Fig. 18 for pure nickel exhibits a negative 
pressure coefficient because the fixed field (/7 
=14) is well above the reversal point, which 
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occurs near H=2. As Yeh has already found? 
the pressure effect in pure iron is negative and 
linear near 71=2. Two other specimens, 20 
percent and 80 percent nickel, also exhibit linear 
changes. See Figs. 9, 10, 15. 

The 30 percent Ni specimen exhibits nearly 
100 percent decrease in flux density at 12,400 
kg/cm? (Fig. 11). This alloy lies near the high 
nickel border of the region in the phase diagram 
of the iron nickel alloys‘ in which the mixed 
crystals of the solid solutions a; and a, can exist 
together. It also lies in the region where the 
Curie point chart® shows a sharp minimum. 
30 percent Ni is about the concentration of the 
alloys which are nonmagnetic, with proper 
preparation, at room temperature. Fig. 3 shows 
that a pressure of 12,400 kg/cm? is sufficient to 
wipe out nearly all the flux density, making the 
metal a very weak ferromagnetic with a per- 
meability about 1.7 which remains practically 
constant from 7=0 to H=100. In Fig. 11 the 
percentage flux change with pressure as the field 
is held fixed at 7.48 is plotted. Unfortunately, 
much trouble was experienced with the pressure 
cylinder during the testing of this specimen, 
preventing a more complete investigation. The 
curves for the first two runs coincide quite well 
and probably indicate the true behavior of the 
metal. Both runs were terminated at 10,500 
kg/cm? by accidental release of pressure. The 
first release was caused by a slow leak in the 
gauge plug packing—the second, by a violent 
explosion of the entire gauge plug. The pressure 
cylinder proved to be defective. Since the state 
of interval stress seems to be closely related to 
the magnetic quality of a ferromagnetic metal, 
an effect of the sudden fall of pressure was 
anticipated. The specimen was accordingly re- 
examined in another press and the third run 
was obtained which indicates an abrupt disap- 
pearance of ferromagnetism near 3200 kg/cm’: 
If this curve represents a true state of the 
metal it must have been induced by the explosive 
release of pressure. In the neighborhood of 30 
percent Ni, the specific volume at 0°C of the 
or nonmagnetic solid solution is less than that 


4 International Critical Tables, Vol. II, p. 451. 
® Invar and Related Nickel Steels, Bur. Standards Circu- 
lar, No. 58, p. 23, Fig. 12. 
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of the magnetic (a:+a2) mixture.’ In_ this 
specimen, therefore, pressure should favor the 
transformation to a nonmagnetic state. 

At low fields 90 and 100 percent Ni reveal 
considerable hysteresis (Figs. 16, 17). That the 
flux does not return to its initial value on release 
of pressure in the case of 90 percent nickel is 
probably due to an inclusion of part of the 
initial irreversible flux decrease. At a higher 
field, H=14, the behavior of nickel reveals a 


® Reference 5, Fig. 41, p. 62. 


curious hysteresis loop (Fig. 18). The mere 
reversal of sign of the pressure increment seems 
to cause a large additional flux change which is 
in the direction in which the flux had just been 
changing. 

In Fig. 20 are plotted the reversible pressure 
changes of flux in percent per 1000 kg/cm? as 
determined by taking differences between the 
“12,000 kg’ and Return to one Atmosphere 
constant pressure curves of Figs. 1 to 8. The 
top curve represents data from Yeh’s paper. 
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CONCLUSIONS 


For all the alloys except 30 percent Ni, the 
following principal effects are to be noted. 

(1) The first application of pressure causes a 
relatively large change, usually a decrease, in 
flux density, followed by incomplete recovery 
on release of pressure. 
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(2) The recovery is reversible in the sense that 
the flux changes give the same curve or hysteresis 
loop regardless of the number of times the 
pressure is applied, provided only that the 
measuring pressures never exceed the first or 
seasoning pressure. 

(3) Reversible changes in flux are in many 
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0 —6 (0k the pressure is changing, except at the terminal 
Tp ‘ie pressures where the two loop sides coincide. 
ye 30” Nickel (4) The rise of the B-JJ curve of a well 
H=0.24 B=926 4 annealed specimen is much sharper than that of 
x First Reversal one not annealed (Fig. 19). (This result, of 
«§ OConstant Reversal — course, is not new.) 
at (5) The effect of a good anneal in changing 
SiS will | = l the flux density is opposite to and greater than 
he Fic. 16. the effect of a 12,000 kg/cm? application of 
he pressure. 
or cases accompanied by hysteresis. Two values of (6) The pressure effect may be positive or 


flux density are found for one given value of negative, depending upon the field strength 
ly . pressure, depending upon the direction in which applied. 
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(7) The field value at which the sign reversal 
in pressure effect occurs depends upon: (a) The 
temperature at which the measurements are 
made. (b) The extent of cold working to which 
the specimen had been subjected. (c) The kind 
of annealing process. 

The results listed above suggest some kind of 
rearrangement of stress pattern within the 
crystal aggregate. It is well known’ that the 
magnetic quality of a ferromagnetic metal is to 
a great extent dependent on the condition of 
internal strain existing within the metal at the 
time of testing and the work of Ruder® indicates 
that it is impossible to remove entirely all the 
strains in a specimen of steel by heating to 
temperatures less than the melting point. If 
such strains do exist in the metal, pressure may 
cause a shift in the strain pattern with accom- 
panying changes in magnetic flux density similar 
to the effects observed by Fraichet® in a magnetic 
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7T. D. Yensen, Magnetic and Other Properties of Iron- 
Aluminum Alloys Melted in Vacuo, Univ. Ill. Eng. Exp. 
Sta. Bull. No. 95, p. 30. 

®W. E. Ruder, Grain Growth in Silicon Steel, Am. 
Inst. Min. & Engrs. Trans. 42, 575 (1913). 

®L. Fraichet, Essai Magnetique des aciers a la traction 
limites elastique, Revue de Metallurgie 20, 549-559 (1923). 
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examination of iron and steel bars while sub- 
jected to simple tension. He found that with 
the first application of tension the flux density 
begins to decrease and falls rapidly and nearly 
linearly with increasing stress with a few super- 
posed small discontinuities until a definite flux 
minimum is reached. Thereafter the flux rapidly 
recovers its original value. The recovery is 
however, not as smooth as the first decrease but 
is accompanied by many small jumps. The flux 
minimum is associated with an elastic limit 
(about half the usual engineering value) which 
Fraichet calls the ‘“‘true elastic limit’? above 
which internal readjustment begins. The flux 
recovery seems to be definitely associated with 
partial release of internal stress by this re- 


adjustment. Fridtjof Schmidt'® has more recently 


1 Fridtjof Schmidt, Investigation of the Induction of 
Iron in Tension Experiments, Diss. Dresden, p. 475 (1931); 
Phys. Bericht. 12, 1247 (1931). 
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carried out similar experiments on magnetic 
changes in steel while under tension. 

Flux measurements such as recorded here are 
the net external effect of an internal magnetic 
field pattern of great complexity. The experi- 
ments of Bitter" on the external magnetic forces 
at the surface of single ferromagnetic crystals 
give information on the fine structure of the 


uF, Bitter, Experiments on the Nature of Ferromagnetism, 
Phys. Rev. 41, 507 (1932). 
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field and show the complexity of the magnetic 
field even in the ideally simple structure of 
single crystals. Because of the complexity of the 
problem a rational explanation of the curves 
submitted would be difficult to formulate. 

The writer wishes to express his indebtedness 
to Dr. L. W. McKeehan for the preparation and 
chemical analyses of the alloys and to Professor 
P. W. Bridgman for full privileges in his labora- 
tory. 
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Direct Measurement of the Gravitational Effect of the Moon 


KENNETH HartLey, Hartley Gravity Balance Corporation, Houston, Texas 
(Received November 30, 1932) 


The sensitivity of the gravity balance described in the 
March, 1932, number of Physics is such that a correction 
must be made for the effect of the attraction of the moon. 
This lunar effect has been measured by hourly readings 
over one complete cycle and agrees, within the observa- 
tional error, with the effect calculated for the vertical 
component of the pull of the moon without any compen- 
sation for the centrifugal effect of the earth’s motion 


around the center of gravity of the system. The effect 
has a period of 24.8 hours and not 12.4 hours as the tides, 
there is no trace of the twelve-hour period. Values of 
gravity plotted against time follow accurately a sine curve 
which shows a small time lag as compared with the 
theoretical curve. This does not seem to be accidental but 
has not been explained. 


N using the gravity balance described in the 
March, 1932, number of Physics for accurate 
measurements, it is necessary to make a cor- 
rection for the effect of the moon. Our unit of 
measure is the millidyne, one-thousandth of a 
dyne per gram or 0.001 cm per sec. acceleration, 
and the effect of the moon when rising or setting 
changes at the rate of about 0.6 millidyne per 
hour, so that in checking back on a station five 
or six hours later than the original observation 
there might be a difference of three millidynes 
or more. I made a rough computation of this 
effect, using round numbers for distance and 
relative mass of the moon and the earth, for 
temporary use, and then decided that it would 
be worth while to make an actual measurement 
of it, keeping the balance in the same place for 
twenty-four hours or more. 

Our location is rather close to the business 
section of Houston, and vibration from street 
traffic is bad, but a basement room under the 
rear end of the building proved much steadier 
than the ground floor so the measurement was 
attempted there. The first intention was merely 
to check the total difference from maximum to 
minimum, but readings in the basement room 
proved so Satisfactory that it was decided to 
try hourly readings and see how closely the 
instrument could follow the gradual change. 

When the ground is fairly steady it is possible 
to read the values of gravity quite definitely to 
one-fifth of a millidyne, but successive readings 
showed considerable irregularity. This was 


proved to be due to imperfect levelling. The 
instrument is intended for field use and so is 
mounted on a tripod and will not remain in 
adjustment very long but must be relevelled for 
each observation. It was found that a variation 
of one-fifth of a millidyne would be produced by 
a disturbance in level of not much more than 
one second, while we cannot level the present 
instrument closer than about five seconds. As the 
maximum variation between several successive 
readings did not exceed the amount attributable 
to imperfect levelling, it is believed that this is 
the only source of error, provided that the 
temperature is constant. In order to eliminate 
this error as far as possible the procedure adopted 
was to take five successive readings, relevelling 
each time, and use the mean of the five as the 
value for the time of the third reading. The 
results by this method were surprisingly close 
even when the deviation of individual readings 
was as much as 1.4 millidynes. 

In the accompanying diagram the solid curve 
is a sine curve with a period of 24.8 hours and 
an amplitude equal to the measured variation 


_in gravity. This was first drawn on transparent 


paper and fitted over the plotted points, as the 
best way to locate the exact time of maximum 
and minimum values, then redrawn on the 
diagram in that position. The points fit the 
curve so closely that it is impossible to shift it 
in either direction more than two or three 
minutes of time without noticeably departing 
from the logical position, yet the times of 
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maximum and minimum do not agree with the 
times of transit of the moon across our meridian 
by about twenty-four minutes. The broken line 
in the diagram is a curve computed from astro- 
nomical data and based on the supposition that 
the instrument is measuring the total pull of the 
moon and not merely the differential effect which 
causes the tides. This differential effect would 
be only a fraction of the observed variation and 
would have a period of twelve hours instead of 
twenty-four. There is no suggestion of the 
twelve-hour period in this curve. 

A preliminary set of readings was made on 
May 6, 1932, not as complete as the later set 


but enough to show the time of maximum and 
minimum quite closely, and the time lag was 
only about three minutes. This suggests that 
the larger lag on May 13 is due to the different 
position of the sun, but the computed effect of 
the sun is much too small to account for it. 

It seems to me that a more extended series of 
observations, perhaps covering a full month or 
more, with a specially built instrument of greater 
precision, might give some very interesting 
results, but I am not in a position to do this so 
thought best to report these preliminary results 
without further delay. 
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The Gliding of a Snake 


ARTHUR TABER JONES, Smith College 
(Received February 15, 1933) 


The snake is replaced by a model consisting of a number 
of spools separated by light pins and held together by 
rubber bands. The forces that act on one spool are con- 
sidered. If the spool is in equilibrium between curved 
guides, and if the tensions in the rubber bands are altered 


in such a way that there is no change in the resultant 
longitudinal forces, but that the curvature of the given 
part of the model tends to increase, it is shown that there 
is a tendency to slip along the guides toward a region 
where the curvature is greater. 


EVERAL recent papers! deal with the gliding 

of snakes. A paper by Fokker, involving 
shearing stresses, is the only one in which the 
dynamical treatment is satisfactory. The fol- 
lowing treatment does not involve shearing 
stresses, and it is simpler than Fokker’s. 

Mosauer’s experimental work shows that when 
a snake is moving with the characteristic gliding 
motion which he calls ‘horizontal undulatory 
movement” the forces which the snake exerts 
on its surroundings are practically normal to its 
body. To simulate the path of a snake among 
grasses, twigs, or stones let a glass tube, large 
enough to admit the snake, be bent into a 
sinuous curve, and let the entire tube lie in a 
horizontal plane. If the snake-glass contact is 
really smooth, can the snake glide through the 
tube? If the snake’s muscles are contracted in 
such a way that he naturally takes the shape of 
the tube he will be in equilibrium in it. If he 
changes his natural curvature to one which 
exists slightly farther along the tube, it seems 
likely that he will slip toward that part of the 
tube. The present note shows that this expecta- 
tion is correct. 

Let the snake be represented by a number of 
spools placed end to end. Let the spools be held 
together by rubber bands, and kept at a fixed 
distance apart by light pins that bear in conical 
depressions in plates on the ends of the spools. 
Fig. 1 shows the forces that act on one spool. 


1A. D. Fokker, Physica 7, 65, 120 (1927). H. J. van 
Leeuwen, Physica 7, 119, 121 (1927). Eduard Wiedemann, 
Zoologische Jahrbiicher, Abt. f. allgem. Zool. u. Physiol. 
d. Tiere 50, 557 (1932). Walter Mosauer, ibid. 52, 191 
(1932); Science 76, 583 (1932). 


Fic. 1. The forces that act on one spool. 


The 7’s are exerted by the rubber bands, the: 
P’s by the pins, and the F’s by the glass tube. 
If the spools are almost in contact on the inner 
side of the curve the angles marked with the 
same symbols will be nearly equal. 

Take the origin at O. Let the mass of the 
spool be m, and its radius of gyration about the 
Z axis be k. Write the three equations which 
state that the sum of the X components of the 
forces equals m#, the sum of the Y components 
equals mj, and the sum of the torques about the 
Z axis equals mk?s—the positive sign for the 
torques being that indicated in the figure. Let 
(1) stand for these three equations when the 
spool is in equilibrium. In this case #= j= d= 0. 

Now let the tensions in the rubber bands be 
changed in such a way that the resultant of the 
forces which the adjacent spools exert on the 
given spool is not changed, but that the torques 
which they exert on it are changed. To accom- 
plish this let 7, and 7: become respectively 
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T: +e and and become 
and 72’—e, P; and P2 be unchanged, and F, and 
F, become F,+e,’ and F2+e’. Let (2) stand 
for the three equations that correspond to (1) 
under the new conditions. If the spool is no 
longer in equilibrium it will start to move, and 
at this instant we shall have j=0, 
Subtract the three Eqs. (1) from the three 
Eqs. (2), and from the three resulting equations 
eliminate ¢,’ and €2’. We obtain 


mé{a+(26 cos 6; cos 92)/(sin (0:+62)) } 
= mk?w+ 2ea(cos 6;—cos 62). (3) 
This equation involves both # and o. To 
obtain another relation between these quantities 
consider the rotation of the given spool when it 
has moved to the right far enough to take the 
place of the next spool. If c stands for the length 
of a pin, and ¢ for the angle through which the 
spool turns, we have approximately 
x=2b+c, g=—2h, x=}#?, and y= 


From these four equations eliminate x, ¢, and 
t, and we see that 


—[262/(2b+c) (4) 


From (3) and (4) eliminate w, and we get 


2b cos 0; cos 62 
sin (0;+62) 
= 2ea(cos 6;—cos 62). (5) 


If the part of the tube in the neighborhood of 
the given spool either is straight or forms an 
arc of a circle, we have 6.=6,, and therefore 
from (5) #=0. So if the whole tube is straight 
or circular, and if it is really smooth, the snake 
will be unable to work its way through it. The 
snake-glass contact is of course not really 
smooth. In fact, it is stated? that a snake can 
climb a vertical glass wall. But let us continue 
with the assumption that the contact is smooth. 

If 62>, (5) shows that # has the same sign 
as e. That is, if the muscular action of the snake 
tends to produce a greater curvature in any 
short length of its body, it also produces a 
tendency for that part to slip toward a neigh- 
boring region of the tube where the curvature 
is greater, and vice versa. We see also that ¢ is 
proportional to (cos 6;—cos 62). If the curvature 
of the tube is the same at both ends of the 
spool we must have 6.=6;. The more rapidly 
the curvature changes as we pass along the tube 
the greater is the difference between 6. and 4;. 
So the force that urges the snake along the 
tube is greatest where the curvature is changing 
most rapidly. This result agrees with that 
obtained theoretically, but quite differently, by 
Fokker, and it is also in accord with the result 
obtained experimentally by Mosauer. 


2H. F. Hutchinson, Nature 20, 528 (1879). 
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